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The temperature and magnetic field dependence of the in-plane thermal conductivity k of the iron- 
arsenide superconductor Ba(Fei_ a; Co :! ;)2As2 was measured down to T ~ 50 mK and up to H = 15 T 
as a function of Co concentration x in the range 0.048 < x < 0.114. In zero magnetic field, a negligible 
residual linear term in k/T as T — > at all x shows that there are no zero-energy quasiparticles 
and hence the superconducting gap has no nodes in the ab-plane anywhere in the phase diagram. 
However, the field dependence of k reveals a systematic evolution of the superconducting gap with 
doping x, from large everywhere on the Fermi surface in the underdoped regime, as evidenced by 
a flat k(H) at T — > 0, to strongly fc-dependent in the overdoped regime, where a small magnetic 
field can induce a large residual linear term, indicative of a deep minimum in the gap magnitude 
somewhere on the Fermi surface. This shows that the superconducting gap structure has a strongly 
fc-dependent amplitude around the Fermi surface only outside the antiferromagnetic/orthorhombic 
phase. 

PACS numbers: 74.25. Fy, 74.20.Rp,74.70.Dd 



The pairing symmetry of the order parameter in iron- 
arsenide superconductors remains a subject of debate as 
different experimental probes point to different conclu- 
sions about the gap structure [l[. For example, while 
the power-law dependence of the penetration depth Q 
and rapid increase of thermal conductivity with magnetic 
field H in (Bai_ x K E )Fe 2 As2 (K-Bal22) can only be ex- 
plained if the gap is very small on some portion of the 
Fermi surface, ARPES studies on the same material find 
a large gap amplitude everywhere [H, [1| , in the experi- 
mentally accessible portions of the Fermi surface close to 
k z =0. Some of the apparent contradictions may come 
from the fact that most studies so far have been carried 
out in a region of the phase diagram close to the anti- 
ferromagnetic (AFM) phase, introducing an additional 
level of complexity. One of the few iron-arsenide systems 
which at present allows a systematic study over the whole 
doping range, including the strongly overdoped regime, 
is Ba(Fei_ K Co K ) 2 As 2 (Co-Bal22) 0,0. 

Here we report measurements of the thermal conduc- 
tivity k of Co-Bal22 throughout the superconducting re- 
gion. By going down to 50 mK we are able to detect 
any zero-energy quasiparticles that would exist if there 
are nodes in the gap. We find a negligible residual linear 
term kq/T in k/T as T — > for all doping levels stud- 
ied, which is strong evidence that there are no nodes 
in the superconducting gap of Co-Bal22, for k direc- 
tions in the a6-plane. However, the ease with which a 
magnetic field H applied along the c-axis induces a fi- 
nite kq/T varies dramatically across the phase diagram. 
Our results show that the superconducting gap evolves 



from being uniformly large everywhere on the Fermi sur- 
face at low doping to having a very small value some- 
where on the Fermi surface at high doping. This reveals 
a nodcless yet highly fc-depcndcnt in-plane gap outside 
the AFM/orthorhombic phase. 

Single crystals of Co-Bal22 were grown from 
FeAs:CoAs flux, as described elsewhere [7[. The dop- 
ing level in the crystals was determined by wavelength 
dispersive electron probe microanalysis, which gave a Co 
concentration x roughly 0.7 times the flux load composi- 
tion (or nominal content). We studied four compositions: 
x = 0.048 (underdoped), 0.074 (optimally doped), 0.108 
(overdoped), and 0.114 (overdoped). The superconduct- 
ing transition temperature (where p = 0) for each sample 
is T c = 17.1, 22.2, 14.6, and 10.1 K, respectively. Sam- 
ples were cleaved into rectangular bars with typical size 
2 - 3 x 0.3 x 0.05 mm 3 . Silver wires were attached 
to the samples with a silver-based alloy, producing con- 
tacts with ultra-low resistance (< 100 Thermal con- 
ductivity was measured along the [100] direction in the 
orthorhombic/tetragonal basal (ab) plane in a standard 
"one-heater-two-thermometer" technique. The magnetic 
field H was applied along the [001] direction (c-axis). 
All measurements were done on warming after cooling in 
constant H from above T c to ensure a homogeneous field 
distribution. 

In Fig. [TJ we show the temperature-dependent electri- 
cal resistivity p(T) of our samples. In the x = 0.048 
sample, the upturn below 60 K signals the onset of 
tetragonal-to-orthorhombic structural transition at Ts 
followed by the AFM order at a temperature Tm 0- 
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FIG. 1: Temperature dependence of the in-plane resistivity 
p(T) of the four single crystals of Ba(Fei_ a; Co :E )2As2 stud- 
ied in this work, with Co concentrations x as indicated. The 
curves are normalized to their room-temperature value prt 
and rigidly offset by 0.05 for clarity. Inset: phase diagram 
showing the orthorhombic phase below Tg, the antiferromag- 
netic phase below Tm and the superconducting phase below 
T c (from 0,11), with arrows indicating the composition of the 
four samples. 



No such feature is seen in the other samples. A smooth 
extrapolation of p(T) to T — yields the residual re- 
sistivity po, equal to 203, 62, 59 and 59 /zfi cm for 
x = 0.048, 0.074, 0.108 and 0.114, respectively. p a is 
used to determine the normal-state thermal conductiv- 
ity kn/T in the T — > limit via the Wiedemann- Franz 
law, k n /T = Lo/po, where L = 2.45 x 10" 8 WO/ K 2 . 
Because the same contacts are used for electrical and 
thermal measurements, there is no relative uncertainty 
between the measured k and this electrically-determined 
kn, other than the small uncertainty in extrapolating 
p(T) to T = 0. 

The thermal conductivity k(T) of the undcrdopcd, op- 
timally doped and most overdopcd samples is shown in 
the left panels of Fig. [21 for a magnetic field ranging from 
H = to 15 T. The data below 0.25 K is well described 
by the function k/T = a + bT a , with a w 1 - 1.5. The 
first term, a = KnlT, is the electronic residual linear term 
of interest here Q. The second term is due to phonons, 
scattered mostly by crystal boundaries at those low tem- 
peratures [ll[. We attribute the fact that this second 
term is much smaller in the x = 0.048 sample, even 
though all samples have comparable size, to the pres- 
ence of structural domains produced by the tetragonal 
to orthorhombic transition at this doping [l2j]. The size 
of domains decreases to ~ l close to a compositional 
boundary a; s =0.073±0.03 (at which T s (x) -> 0) [13], and 
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FIG. 2: Left panels: Temperature dependence of the thermal 
conductivity, plotted as k/T vs T, for three Co concentra- 
tions, each measured in magnetic fields as indicated. The 
lines are a power law fit of the form k/T = a + bT a to the 
H = and H — 15 T data, used to extract the residual linear 
term kq/T = a in the T — > limit (see text). Right panels: 
field dependence of k/T plotted for three temperatures, as in- 
dicated. The T — > data (red squares) are obtained from the 
power-law extrapolations, with a typical error bar as shown. 



diminishes effective phonon mean free path. The opti- 
mally doped sample x=0.074 is located at the border of 
the orthorhombic phase and its phonon conductivity is 
still notably suppressed. 

The magnitude of the residual linear terms extracted 
from the fits to the H = data in the left panels of Fig. [2] 
are extremely small, below 5 pW / K 2 cm in all cases. 
This is within the absolute error bar of ± 5 pW / K 2 cm 
on kq/T obtained in similar measurements on samples 
where no residual linear term is expected (llj . Within 
those error bars, all Co-Bal22 samples exhibit a negligi- 
ble residual linear term, irrespective of the doping level. 
Similarly negligible values of kq/T were observed in K- 
Bal22 [3j and, as we argued there, this allows us to safely 
conclude that the superconducting gap of Co-Bal22 docs 
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not contain a line of nodes anywhere on the Fermi sur- 
face, at any doping. By comparison, in a d-wave super- 
conductor like the cuprate Tl 2 Ba 2 CuOfi-i_,i; (Tl-2201) with 
T c ~ 15 K, k /T = 35 % of k n /T while in all four 
Co-Bal22 samples, k q /T < 1 - 2 % of k n /T. This rules 
out d-wave symmetry even for the overdoped samples of 
Co-Bal22 and further confirms that cuprates and iron 
arsenides have distinct superconducting states. 

The possibility of point nodes in the gap of Co-Bal22 
was raised in the context of penetration depth measure- 
ments which gave AA oc T 13 with (3 rj 2 14|. Zero- 
energy quasiparticles associated with symmetry-imposed 
point nodes give a residual linear term in n which grows 
with impurity scattering (if the gap grows linearly from 
the node) and for high scattering rates kq/T should be- 
come a substantial fraction of kn/T [15j. (If the gap 
grows quadratically from the node, kq/T is independent 
of scattering rate and comparable to the case of a line 
node [H}.) The normal-state impurity scattering rate Tq 
of Co-Bal22 can be estimated from po and the plasma 
frequency uj p = c/Ao via Tq = (ojp/87r)/9o, where Ao is 
the penetration depth at T — > and c is the speed of 
light. At optimal doping, Ao — 200 nm 
hYo/k^Tc ~ 1 — 2 in the samples x 
and 0.114. This is very substantial, and would give a 
large kq/T from point nodes, comparable to the case of a 
line node. We therefore conclude that symmetry-imposed 
point nodes in the gap of Co-Bal22 are also ruled out, 
unless they are located in the k z direction (along the c- 
axis). 

In a recent study, a large kq/T was reported fl6| in 
a Co-Bal22 sample with T c and residual resistivity ra- 
tio similar to our x = 0.074 sample, amounting to 60 % 
of k-^/T, compared with less than 1-2 % in our sample. 
This stands in contradiction not only with our results 
on more than 10 crystals of Co-Bal22 at various dop- 
ings, but also with previous results on K-Bal22 0] and 
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Ni-doped BaFe 2 As 2 (Ni-Bal22) 
k /T ~ at H = 0. 

The dependence of k on magnetic field H || c is dis- 
played in the right panels of Fig. [2] The purely elec- 
tronic component is given by the extrapolations to T = 
(T —> 0; red squares). To confirm that the H dependence 
is not affected by the extrapolation procedure, we com- 
pare the extrapolated data with cuts through the actual 
data at two fixed temperatures, namely 60 mK and 100 
mK. We see that the latter cuts produce what is essen- 
tially a rigid shift of the data relative to T — » 0, caused 
by the added phonon contribution. With this validation 
of the T — > data, we plot it on normalized scales in 
Fig.[3]as (k /T)/(k n /T) vs H/H c2 , where H c2 is the up- 
per critical field, obtained from resistivity measurements 
in high fields @, Oj, which give H c2 ~ 45, 60, 30, and 
20 T for x = 0.048, 0.074, 0.108, and 0.114, respectively. 
This plot reveals the ease with which zero-energy quasi- 
particles are excited by a magnetic field. At x = 0.048, 




FIG. 3: Left panel: Residual linear term ko/T in the thermal 
conductivity of Co-Bal22 (full symbols) as a function of mag- 
netic field H, plotted on normalized scales as {ko/T)/{kn/T) 
vs H/H C 2, where kn is the normal-state thermal conductiv- 
ity in the T — > limit, obtained from the Wiedemann- Franz 
law (see text), and H C 2 is the upper critical field, for the four 
Co concentrations x as indicated. For comparison, we also 
show the behaviour of pure (Nb) and dirty (InBi) s-wave su- 
perconductors (from [3]), and of the rf-wave superconductor 
Tl-2201 {T c = 15 K and H c2 ~ 7 T) (from pjjj). The data 
for the borocarbide superconductor LuNi2B2C (T c = 16 K 
and H C 2 = 7 T) is also shown (from [2TJ]). Right panel: 
(reo/T) /(kn/T) at H/H c2 = 0.25 in Co-Bal22 as a function of 
x. This measures the extent to which zero-energy quasiparti- 
cles are excited by a relatively small magnetic field compared 
to H C 2- The Co-Bal22 data shows a degree of excitation 
typical of s-wave superconductors in the region of coexisting 
antiferromagnetic/orthorhombic order at low doping, rising 
rapidly to a level typical of a d-wave superconductor at higher 
doping. 



Co-Bal22 shows the slow rise with upward curvature typ- 
ical of isotropic s-wave superconductors like Nb or InBi 
(see Fig. [3]), where the rise in kq/T is exponentially slow 
because it relies on the tunneling of quasiparticles be- 
tween localized states inside adjacent vortex cores, which 
at low fields are far apart. By contrast, at x = 0.114, Co- 
Bal22 shows a very rapid rise with downward curvature 
more typical of d-wave superconductors such as Tl-2201 
(see Fig. Even though the superconducting gap does 
not go to zero anywhere, as established by the fact that 
kq/T ~ 0, it must become very small somewhere on the 
Fermi surface to account for such a rapid rise with field. 
We can estimate roughly the minimum value of the gap 
A from the fact that at x = 0.114 the rise in Ko/T from 
zero starts around H ~ 0.5 T (see lower right panel in 
Fig.HJ), i.e. at H/H c2 ~ 1/40. Given that H c2 oc A 2 , we 
obtain for the gap minimum A m ; n ~ Ao/6, in terms of 
the gap maximum, Aq, which controls H c2 . 
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This is similar to what has been found in the borocar- 
bide superconductor LuNi 2 B 2 C (T c = 16 K and H c2 = 
7 T), where even though kq/T = at H = 0, k /T rises 
very rapidly at low fields [20j (see Fig. [3]). There, an 
estimate of the gap minimum gave A m ; n ~ Ao/10 [20|. 
Recent ARPES studies on the closely related material 
YNi2B2C found a sharp superconducting gap minimum 
at positions on the Fermi surface which are connected by 
a nesting wavevector 
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The fact that the deep minimum in the gap of Co- 
Bal22 is only present in composition free of competing 
orders (in our case orthorhombic/AFM phase) is consis- 
tent with previous heat transport studies performed at 
single dopings on other arsenide superconductors. There 
is no AFM order in K-Bal22 at x ~ 0.3, where a rela- 
tively rapid increase of kq/T vs H was reported in- 
dicative of a gap minimum. In the low-T c nickel-arsenide 
superconductor BaNi 2 As 2 (T c = 0.7 K) [23 |. which or- 
ders around 132 K with an as yet unknown order pa- 
rameter [23| , the weak field dependence of k is consistent 
with a full gap everywhere on the Fermi surface, com- 
patible with the behavior in undcrdoped Co-Bal22. In 
Ba(Fci_ x Ni :c )2As2 (Ni-Bal22), only a weak H depen- 
dence was observed near optimal doping (T c = 20 K) 
fl7| - it will be interesting to see how that evolves at 
higher doping. Note that in Ni-Bal22, the competing 
magnetic /structural phase may extend to higher electron 
doping [241 ] . In the low-T c iron-arsenide superconductor 
LaFePO (T c = 7 K), a line of nodes in the gap has been 
inferred from the linear T dependence of the penetration 
depth Such a line of nodes would necessarily lead 

to a finite kq/T, and indeed it does, as reported recently 
[2^ . This extreme form of anisotropy again occurs in the 
absence of competing orders. 

In summary, zero-field heat transport in the T — > 
limit shows that the superconducting gap of Co-Bal22 
does not go to zero at any point on the Fermi surface, 
for quasiparticle k vectors in the basal plane. This is 
true throughout the phase diagram, from underdoped to 
overdoped, and is consistent with heat transport data 
in both K-Bal22 and Ni-Bal22 [Hi, measured previ- 
ously at a single doping close to optimal. In Co-Bal22, 
the initial rise of kq /T with magnetic field H varies enor- 
mously with doping, from slow and typical of an isotropic 
s-wave gap at low doping, in the region of coexisting 
AFM/orthorhombic order, to very fast and comparable 
to a d-wave gap at high doping. This fast initial rise 
shows that the gap must be very small on some portion 
of the Fermi surface, roughly a sixth of the gap maximum. 
One way to understand this pronounced fc-dependence of 
a nodeless super conduc ting gap is in terms of an extended 
whereby the symmetry of the or- 
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gap zeros or nodes are 'accidental', i.e. not imposed by 
symmetry (as they would be for a d-wave order param- 
eter) , the effect of impurity scattering is not to broaden 
them (as in a d-wave gap), but to lift them, making the 
gap more isotropic [28| . This could then explain why 
the undoped material LaFePO has a line of nodes, while 
the doped superconductor Co-Bal22 doesn't have nodes, 
given that the former has a 10 times smaller impurity 
scattering rate (as measured by the residual resistivity 
ratio). In other words, as recently proposed j^, the pro- 
nounced fc-dcpcndcncc of the gap produces nodes in the 
pure limit which turn into deep minima in the presence 
of strong impurity /disorder scattering. The fact that no 
pronounced anisotropy (i.e. field dependence of re) is ob- 
served in Co-Bal22 when there is AFM order leads us 
to propose that the AFM order in iron arsenides recon- 
structs the Fermi surface by gapping those regions where 
the minima in the superconducting gap are located. This 
is reminiscent of what occurs in the borocarbide super- 
conductors i?Ni2B2C (R = Y or Lu) [2(j, where gap min- 
ima are observed at fc-locations connected by a nesting 
wavevector 2l|. This suggests an intimate interplay of 
magnetism and superconductivity 

Work at the Ames Laboratory was supported by the 
Department of Energy-Basic Energy Sciences under Con- 
tract No. DE-AC02-07CH11358. R.P. acknowledges sup- 
port from the Alfred P. Sloan Foundation. L.T. acknowl- 
edges support from CIFAR, NSERC, CFI and FQRNT. 



der parameter is s-wave but the pairing interaction is so 
strongly fc-dependent that the gap is highly anisotropic, 
going to zero in certain directions. (In addition, it may 
also change sign, as in the s± state |30|.) Because these 



[2 
[3; 

[5 
[6 
[7 
[8 
[9 
[10 

I" 
[12 
[13 

[14 
[15 
[16 

[17 

[18 
[19 
[20 
[21 
[22 
[23 
[24 



E-mail: louis.taiUefer@physique.usherbrooke.ca 

K. Ishida, Y. Nakai, and H. Hosono, J. Phys. Soc. Jpn. 

78, 062001 (2009). 

C. Martin et al, Phys. Rev. B 80, 020501(R) (2009). 
X.G. Luo et al, Phys. Rev. B 80, 140503 (R) (2009). 
H. Ding et al, Europhys. Lett. 83, 47001 (2008); L. Zhao 
et al, Chin. Phys. Lett. 25, 4402 (2008). 

K. Nakayama et al, Europhys. Lett. 85, 67002 (2009). 
A.S. Sefat et al, Phys. Rev. Lett. 101, 117004 (2008). 
N. Ni et al, Phys. Rev. B 78, 214515 (2008). 

D. K. Pratt et al, Phys. Rev. Lett. 103, 087001 (2009). 
H. Shakeripour et al, New J. Phys. 11, 055065 (2009). 
R. Prozorov et al, Phys. Rev. B bf 80, 174517 (2009). 
S.Y. Li et al, Phys. Rev. B 77, 134501 (2008). 

M.A. Tanatar et al, Phys. Rev. B 79, 180508 (R) (2009). 

C. Proust et al, Phys. Rev. Lett. 89, 147003 (2002). 

R.T. Gordon et al, Phys. Rev. Lett. 102, 127004 (2009). 

M.J. Graf et al, Phys. Rev. B 53, 15147 (1996). 

Y. Machida et al, J. Phys. Soc. Jpn. 78 073705 (2009). 

L. Ding et al, New J. Phys. 11, 093018 (2009). 

M. Kano et al, J. Phys. Soc. Jpn. 78, 084719 (2009). 

S.Y. Li et al, Phys. Rev. Lett. 99, 107001 (2007). 

E. Boaknin et al, Phys. Rev. Le tt. 87, 237001 (2001). 
Baba et al. ! arXiv:0812.0406l 

Kurita et al, Phys. Rev. Lett. 102, 147004 (2009). 
S. Sefat et al, Phys. Rev. B 79, 094508 (2009). 
C. Canfield et al, Phys. Rev. B 80, 060501 (2009). 



5 



[25] J.D. Fletcher et ai, Phys. Rev. Lett. 102, 147001 (2009). 
[26] M. Yamashita et ai. larXiv:0906.0622l 
[27] F. Wang et ai, Phys. Rev. Lett. 102, 047005 (2009). 
[28] V. Mishra, et ai, Phys. Rev. B 79, 094512 (2009). 



[29] A.V. Chubukov et ai, Phys. Rev. B 80, 140515(R) 
(2009). 

[30] I.I. Mazin et ai, Phys. Rev. Lett. 101, 057003 (2008). 



